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In the presence of strong magnetic field reported to have been observed on the 
surface of some neutron stars and on what are called Magnetars, a host of physical 
phenomenon from the birth of a neutron star to free streaming neutrino cooling 
phase will be modified. In this review I will discuss the effect of magnetic field on 
the equation of state of high density nuclear matter by including the anomalous 
magnetic moment of the nucleons into consideration. I would then go over to dis- 
cuss the neutrino interaction processes in strong as well as in weak magnetic fields. 
The neutrino processes are important in studying the propagation of neutrinos and 
in studying the energy loss, Their study is a prerequisite for the understanding of 
actual dynamics of supernova explosion and on the stabilization of radial pulsation 
modes through the effect on bulk viscosity. The anisotropy introduced in the neu- 
trino emission and through the modification of the shape of the neutrino sphere 
may explain the observed pulsar kicks. 

1 Introduction 

Large magnetic fields ~ 10 14 Gauss have been reported to exist on the sur- 
face of pulsars. Recent observations ID of 7-ray repeaters and spinning X-ray 
pulsars (Magnetars) hints to the existence of fields in excess of 10 14 Gauss. It 
then follows from the Scalar Virial theorem that the fields in the core could 
even reach a value as large as 10 18 Gauss. There is however, an upper limit on 
the magnetic field discussed by Chandershckhar beyond which the magnetic 
energy exceeds the gravitational energy and the star is no longer stable. In the 
presence of such magnetic fields, neutron star properties in all its phases from 
its evolution from proto-neutron star to cold neutrino emitting phase would 
be modified. This arises because in the presence of magnetic field motion of 
the charged particles is quantised in a plane perpendicular to the magnetic 
field and the charged particles occupy discrete Landau levels. This has the 
effect of not only modifying the energy eigenvalues but also the particle wave 
functions. The quantum state of a particle in magnetic field is specified by 
momentum components p z ,p y , spin s and Landau quantum number v. The 
anomalous magnetic moments of protons and neutrons further modify the en- 
ergy eigenvalues. The time scales involved during all phases of neutron stars 
from birth to neutrino burst through thermal neutrino emission from trapped 
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neutrino sphere to freely streaming neutrino cooling phase are large compared 
to the interaction time scales of strong, electromagnetic and weak interactions 
and the matter is in /? equilibrium. The magnetic field would modify the equi- 
librium and all neutrino interaction processes including scattering, absorption 
and production. 

The stategy then is to first solve the Dirac equation for all particles in 
magnetic field including their anamolous magnetic moments, obtain the en- 
ergy eigenvalues, construct the Grand Partition Function taking into account 
strong interactions in some model dependent way and obtain the equation of 
state (EOS). The next step is to calculate the scattering cross-sections for all 
neutrino processes by using the exact wave functions and by modifying the 
phase space integrals for arbitrary values of degenracy, density, temperature 
and magnetic field. The various phenomenon that I will address are : 

1. Composition of matter in neutron stars, proton fraction, effective nucleon 
mass etc. 

2. Cooling of neutron stars in the free streaming regime. 

3. Neutrino transport in neutron and collapsing stars which is an essential 
prerequisite for an understanding of supernova explosion, structure of 
proto- neutron star and observed pulsar kicks. 

4. Damping of radial oscillations and secular instability through the calcu- 
lation of bulk viscosity. 



2 Nuclear Matter Composition 



For determining the composition of dense, hot, magnetised matter, we employ 
a relativistic mean field theoretical approach in which the baryons (protons 
and neutrons) interact via the exchange of a — lo — p mesons in a constant uni- 
form magnetic field. Following reference^lin a uniform magnetic field B along 
z axis corresponding to the choice of the gauge field = (0, 0, xB, 0),the rel- 
ativistic mean field Lagrangian can be written as 



i^D 1 * - m B + g<yB<J - 9ujB1^ - 9 p b t 3B 7m P M + kbct^F^ 



-m-p^ - U(a) - \f^f^ + - 



4>B 



(1) 
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in the usual notation with = d^ + ieA^ and Kb as the anomalous magnetic 
moments given by 



2m p 



■1L _ i 
2 



2m n 2 



''" (2) 



where g p — 5.58 and g n = —3.82 are the Lande's g-factor for protons and 
neutrons respectively. Replacing the meson fields in the relativistic mean field 
approximation by their density dependent average values < oq >,< uoq > and 
< po >, the equations of motion satisfied by the nucleons in the magnetic 
field become 
d 

[-iotx-x- + a y(Py - eBx ) + a zPz + Pim - 9a < cr >) - iKia x a y B}f itS 

Ox 

= (E i -Ui)f i , s (x) (3) 
where 

K n = g uN <^o> ± 9j f <Po> (4) 

The equations are first solved for the case when momentum along the mag- 
netic field direction is zero and then boosting along that direction till the 
momentum is p z . For the neutrons and protons we thus get 



C = \jmf +p^ + kIB 2 + 2 Kn Bs^ P l+pl+m* 2 

= K - US (5) 



and 



Sis = \jmf+pl + eB(2v + 1 + s) + K 2 p B 2 + 2K v Bs^mf + eB{2v + 1 + s) 

= K, s - u p (6) 

where 2v = (2n+ 1 + s), where v and n being integer known as Landau quan- 
tum number and principal quantum number respectively, s = ±1 indicates 
whether the spin is along or opposite to the direction of the magnetic field 
and m* B = raj — g a a is the effective baryon mass. The energy spectrum for 
electrons is given by 



K,s = y/ml+pl + (2v+l + s)eB (7) 

The mean field values < oo >, < ojq > and < p > are determined by 
minimizing the energy at fixed baryon density rig = n p + n n or by maximizing 
the pressure at fixed baryon chemical potential (is- We thus get 

< a >= -^r«+<) , < >= ^{n n +n p ) , < Po >= --^(n p -n n ) (8) 

777 „ 77%. , A 777 n 
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where rij and nf (i=e,p) are the number and scalar number densities for 
proton and neutron. In the presence of the magnetic field phase space volume 
is replaced as 



dp z 



(9) 



The expressions for number densities, scalar densities for neutron and proton 
are given by 

1 V f d 3 p 

>7T) 3 ^ J 1 



n n = 



n v = 



(2tt) 3 11 ePW-O 

x ' S 

(2tt) 3 ^ J £2(\ + ePW-0) 

eB^^ f°° dp z 



n', 



+ e /3(£?, s -M;) 

dp z 



2n 2 - P ^^j £PvA1 + e m%-^) ) 

and the net electron and neutrino number densities are given by 



n e = 



eB 
2^2 



E( 2 -m / 



dp z 



1 



(1 + e(0(£< - fie))) 



- Me 



(-Me) 



n v . = 



(2tt) 3 



d 3 p 



( 1 + e /J(£"« -*!.,.)) 



M^e <"> (-M^e) 



(10) 

(11) 
(12) 

(13) 

(14) 
(15) 



The thermodynamic potential for the neutron, proton, electron and neutrino 
are given by 



-/3(£ s n -iOl 



^^/*.Ml + e-«M] 



Hp — 



n 



-/J(^,.-Me)l 



+ e' 



-0(£"«-M„ e )l 



(16) 

(17) 
(18) 
(19) 
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and the thermodynamic potential for the system is given by 

n = -\™L^l - \my + \mlo 2 + U(a) + |- + ]T fi 4 (20) 

where i = n,p,e,v e . The various chemical potential are determined by the 
conditions of charge neutrality and chemical equilibrium. In later stages of 
core collapse and during the early stages of protoneutron star, neutrinos are 
trapped and the chemical potentials satisfy the relation \x n = p p + p e — fi Ve 
These situations are characterized by a trapped lepton fraction Yl = Y e + Y Ve 
where Y e = (n e —n e+ )/nB is the net electron fraction and Y Ve — {n Ve —nv e )/nB 
is the net neutrino fraction. 

The evolution of a protoneutron star begins from a neutrino-trapped sit- 
uation with Yl ~ 0.4 to one in which the net neutrino fraction vanishes and 
chemical equilibrium without neutrinos is established. In this case the chem- 
ical equilibrium is modified by setting fi Ve = 0. In all cases, the conditions of 
charge neutrality requires 

n p = n e (21) 

In the nucleon sector, the constants g a N, g^N, gpN, b and c are determined by 
nuclear matter equilibrium density no = 0.16/m~ 3 and the binding energy per 
nucleon (<~ WMeV), the symmetry energy (~ 30 — 35MeV), the compression 
modulous (200MeV < ko < SOOMeV) and the nucleon Dirac effective mass 
M* = (0.6 — 0.7) x 939MeV at n . Numerical values of the coupling constants 
so chosen are : 

^ = 3.434/m- 1 : ^ = 2.694/m- 1 ; ^ = 2.1 fm' 1 
m<j m u m p 

b = 0.00295, c= -0.00107. 
3 Weak Rates and Neutrino Emissivity 

The dominant mode of energy loss in neutron stars is through neutrino emis- 
sion. The important neutrino emission processes leading to neutron star cool- 
ing are the so called URCA processes 

n — >p + e~+D e p + e~H ► n + v e (22) 

At low temperatures for degenerate nuclear matter, the direct URCA process 
can take place only near the fcrmi energies of participating particles and 
simultaneous conservation of energy and momentum require the inequality 

p F (e) +p F (p) > Pf(ii) 
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to be satisfied in the absence of the magnetic field. This leads to the well 
known threshold El for the proton fraction Y p — ^ > 11% thus lead- 
ing to strong suppression in nuclear matter. This condition is satisfied for 
jib > 1.5no ( where no — 0.16/m _3 ,is the nuclear saturation density), in a 
relativistic mean field model of interacting n-p-e gas for B = 0. The standard 
model of the long term cooling is the modified URCA: 

(n,p) + n — > (n,p) + p + e~ + D e , (n,p)+p + e~ — >{n 1 p) + n + v e (23) 

which differ from the direct URCA reactions by the presence in the initial 
and final states of a bystander particle whose sole purpose is to make pos- 
sible conservation of momentum for particles close to the Fermi surfaces. For 
weak magnetic field the matrix element for the process remains essentially 
unaffected and the modification comes mainly from the phase space factor. 
Treating the nucleons non-relativistically and electrons ultra-relativistically, 
the matrix element squared and summed over spins is given by 

\ M ? = 8G f cos 2 c (4m*m;)E e E v [(l + 3g 2 A ) + (1 - g\) cosfl c ] (24) 

where gA = 1.261 is the axial- vector coupling constant. The emissivity ex- 
pression is given by 



£ v = 



E v V \M\\2nf8\P f - Pi)S (25) 



where the phase space integrals are to be evaluate over all particle states. 
The statistical distribution function S = f n (l — f p )(l — f e ), where / z 's are the 
Fermi-Dirac distributions. We can now evaluate the emissivity in the limit 
of extreme degeneracy, a situation appropriate in neutron star cores by using 
the standard techniques to perform the phase space integrals^ 

4577r i 
t = J^° 2 f COs2 W + ^ 2 A )m* n m;eBT e ^ (2 - 5 U . ) == 

(26) 

( '-mA 

where u max = Int Me 2e g' e . In the limit of vanishing magnetic field, the sum 

can be replaced by an integral and we recover usual expression i.e. the case 
for B=0. 

4577T 

£»{B = 0) = i^G 2 F cos 2 O (1 + 3g A )m* n m;p e T 6 (27) 

Modified URCA process are considered to be the dominant process for neutron 
star cooling. Similarly the energy loss expression with appropriate electron 
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phase space, for modified URCA process Surca is calculated to be 

• 11513 G\ cos 2 9 C a .3 / V tB 

Wa= 60480 2n ^"H^J 

'(2-^,0) {-==4== ) (28) 



E 

!/=0 



^/ — to 2 — 2i/ei3 



where / is the ir ~ N coupling constant (/ 2 ~ 1) and Qurca has been esti- 
mated to be = 1.54. The above equation in the B — > limit goes over to the 
standard result 

• 11513 g| cos 2 6 C a » 3 / / V ^8 , 9cn 

WA(B = 0) = 3024Q 2?r gA"*n ™P J MeT (29) 

In the case of super strong magnetic fields such that B > B e c (B e c — 4.41 x 10 13 
Gauss) all electrons occupy the Landau ground state at T=0 which corre- 
sponds to v = state with electron spins pointing in the direction oppo- 
site to the magnetic field. Charge neutrality now forces the degenerate non- 
relativistic protons also to occupy the lowest Landau level with proton spins 
pointing in the direction of the field. In this situation we can no longer con- 
sider the matrix elements to be unchanged and they should be evaluated using 
the exact solutions of Dirac equation. Further because nucleons have anoma- 
lous magnetic moment, matrix elements need to be evaluated for specific spin 
states separately. The electron in v — state has energy E e = \Jm 2 e + p^ z 
and the positive energy spinor in v = state is given by 

( \ 

E e + TO, 



o 

\ -Pez 

Protons are treated non-relativistically and the energy in v — state is E p 
m p + ^4 — h Uq and the non-relativistic spin up operator J7p,+i given by 



(30) 



1 



\/2to p 



= -7= ( X n + ) (31) 







where rh p = to* — k p B. For neutrons we have 



,,,(r) = e- ip ""- r U n ,s(E n ,s) (32) 

Y l-l X LyL Z 

U ns = ^=( X A (33) 
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and E n ^ s ~ m* + j5^ 2 — n n Bs + Uq in the non-relativistic limit. The neutrino 
wave function is given by 

*„(r) = 1 e-'f-U^Ev) (34) 
y L x Ly-L z 

Here Z7„ )S is the usual free particle spinor, %s is the spin spinor and the wave 
function has been normalised in a volume V = L x L y L z . Using the explicit 
spinors given above we can now calculate the matrix element squared and 
summed over neutrino states to get 

Y, \MlM+\ = \K + \H(4m* n m p )(l+g A ) 2 (E e + p ez ){E v+Pvz ) (35) 

and 

^|M!M_| = \K^\ 2 16(Am* n m p )g 2 A (E e +p ez )(E u -p uz ) (36) 

The neutrino emissivity is calculated by using the standard techniques for 
degenerate matter and we get 



S= cos 2 6 c eB^Tt 

40320 F p F (e) 



2 V J \ 2eB 



- [ P%(n, +) - 4 P 2 (e) ) e J zmh3_^m 



+2g 2 APF(n,+) — >p F (n,- 



(37) 



where Pf(n, +) and p F {n 1 —) are the neutron Fermi momenta for spins along 
and opposite to the magnetic field direction respectively and are given by 

-^*r = ^- m »- Uo± ^K (38) 

We thus see as advertised that in the presence of quantising magnetic field 
the inequality p F (e) + p F (p) > Pf(n) is no longer required to be satisfied for 
the process to proceed, regardless of the value of proton fraction and we get 
non zero energy loss rate. 

4 Bulk Viscosity of Magnetised Neutron Star Matter 

The source of bulk viscosity of neutron star matter is the deviation from f3 
equilibrium, and the ensuing noncquilibrium reactions, implied by the com- 
pression and rarefaction of the matter in the pulsating neutron star. These 
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important reactions are the URCA and the modified URCA processes. Since 
the source of bulk viscosity is the deviation from (3 equilibrium, these reactions 
are driven by non zero value of A/i = ji n — ji p — /i e . We calculate the bulk 
viscosity of neutron star matter in the presence of magnetic field for direct 
URCA processes in the linear regime i.e. A/i << kT . The bulk viscosity £ 
is defined bj0 

f -'(S)-=(£)> . <39> 

Here vo is the specific volume of the star in equilibrium configuration, Av is 
the amplitude of the periodic perturbation with period r = — and v(t) = 

v o + Av sin (^—f^j The quantity (^f-)aw is the mean dissipation rate of energy 
per unit mass and is given by the equation 



. , i = — / P(t) — dt (40) 
\dtj av rj V, dt 

The pressure P(t) can be expressed near its equilibrium value Po, as 

p «= p » + (f ) a « + (So an * + (£)o a "- + (S„ a "' (4l) 

The change in the number of neutrons, protons and electrons per unit mass 
over a time interval (0, t) due to URCA reactions (23) is given by 

- An n = An p = An e = [ ^-dt (42) 

The net rate of production of protons, ^f- is given by the difference between 
the rates 7\ and r 2 of the URCA reactions. At equilibrium the two rates 
are obviously equal and the chemical potentials satisfy the equality A/i = 
Mn — fJ-p ~ A*e = A small volume perturbation brings about a small change 
in the chemical potentials and the above inequality is no longer satisfied, now 
Afi is not zero and consequently the reaction rates are no longer equal. The 
net rate of production of protons will thus depend upon the value of Afi. In 
the linear approximation, ^ -C 1, the net rate can be written as 

dr ^ = r 1 -r 2 = -XAn (43) 



dt 

Using the thermodynamic relation — — j^jj- and employing the above 
relation we obtain 

SP = P(t) -P Q = J XAfi(t)dt (44) 
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The change in the chemical potential A/i(t) arises due to a change in the 
specific volume Av and changes in the concentrations of various species, viz, 
neutrons, protons and electrons. Thus 



Afj,(t) = A/i(0) 



1 > Av + — — An, 



V /o 
fdA^\ 



V <9n„ / 



/ 9A/A 



and we arrive at the following equations for Afj,: 



dAfi 
dt 



Av 

uA — cos(wi) — CXAfi 



(45) 



(46) 



where 



4 



/ dA^ 9A^ <9A M \ 
V ° V <9n p <9n e <9n„ / 



(47) 



Since for small perturbations, A, A and C are constants, equation (46) can be 
solved analytically to give 

"~ A A '' C\e- cxt +Lusm(ujt) + C\cos{ujt) (48) 



A^i = 



C 2 \ 2 



and we obtain the following expressions for £ 



c = 



A 2 X 



C 2 X 2 



1 - 



;CA 1 - e 



-CXt 



7T CJ 2 + C 2 A 2 



(49) 



Given the number densities of the these particle species in terms of their 
respective chemical potentials, one can determine the coefficients A and C; 
given the rates Ti and for the two URCA processes one can determine 
A and hence £ for any given baryon density and temperature. For weak 
magnetic field several Landau levels arc populated and the matrix elements 
remain essentially unchanged and one needs to account for the correct phase 
space factor. For non-relativistic degenerate nucleons the decay rate constant 
A is given by^ 



A = 



17 

480tt 

^r 2 



G 2 F cos 2 Bc T 4 (l + 3g%)eB m* p m* n 9 (p F {p) + p F {e) - p F (r 



x 



v=0 



S v ,o] 



IveB 



(50) 
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For strong magnetic field, the electrons are forced into the lowest Landau 
level. Using the exact wave functions for protons and electrons in the lowest 
Landau level and carrying out the energy integrals for degenerate matter, the 
decay constant A is given by© 



1 7 77Z* 771* 

A = ^T G F cos l eB^-^T 4 [Ag 2 A + (g v + g A ) 2 ] 
9607T Pf(P) 



exp 



2eB 



^(/>J. («> - 4y»7,(r)) <r P [ >^± ) t , 



It is clear from above that in the case of completely polarised electrons and 
protons the direct URCA decay rate always gets a non-zero contribution from 
the second term in the last square bracket irrespective of whether the trian- 
gular inequality pf(s) + pf(p) > Pf(it>) is satisfied or not. 

5 Neutrino Opacity in Magnetised Hot and dense Nuclear 
Matter 

We calculate neutrino opacity for magnetised, interacting dense nuclear mat- 
ter for the following limiting cases: a) nucleons and electrons, highly degener- 
ate with or without trapped neutrinos, b) non-degenerate nucleons, degenerate 
electrons and no trapped neutrinos and finally, c) when all particles are non- 
degenerate. The important neutrino interaction processes which contribute 
to opacity are the neutrino absorption process 

v e + n — ► p + e (51) 

and the scattering processes 

v e + N -> v e + N 
v e + e — > v e + e (52) 

both of which get contributions from charged as well as neutral current weak 
interactions. For the general process 

v( P i)+A(p 2 )^B(p 3 ) + l{p i ) (53) 

The cross-section per unit volume of matter or the inverse mean free path is 
given by 



V X ~ 1[El) =2E n d Pi W fi M**) [1-/3(^3)] [1-/4(^4)] 

1 i=2,3,4 



(54) 
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where dpi = ^ 2 t)^2E- lSl ^ ne density of states of particles and the transition 
rate W fi = (2tt) 4 ( 5 4 (P / - P % ) \M\ 2 . 

Weak Magnetic Field : For weak magnetic fields, several Landau levels are 
populated and the matrix element remain essentially unchanged and one needs 
to only account for the correct phase space factor. We first consider the 
neutrino-nucleon processes. In the presence of weak magnetic fields, the ma- 
trix element squared and summed over initial and final spins in the approxi- 
mation of treating nucleons non-relativistically and leptons relativistically is 
given by 



J2\M\ 2 = 32G 2 F cosl m>* 



(<# + SC%) + ((% - C 2 A )cosl 



EpE,, 



(55) 



where CV = gv =1, Ca = 9a =1-23 for the absorption process; Cy — 
— 1, Ca = —1.23 for neutrino scattering on neutrons and CV = — 1 + 
Asin 2 w =0.08, Ca =1.23 for neutrino proton scattering. We now obtain neu- 
trino cross-sections in the limits of extreme degeneracy or for non-degenerate 
matter. 

Degenerate Matter : The absorption cross-section for highly degenerate mat- 
ter can be calculated by using (79) in (78) by the usual techniques and we get 
for small B El 



V 



G F COSg c 2 

8^ 



3g 2 A )m*m* n T 2 



eB 



Pf(p) + PF(e)p F (n) - pf{v) 



[pf{p) +PF{e) -PF{n) +Pf(v)] 
2E U 



Pf{v) - \pf(p) +PF(e) -p F (n)\ 



v=0 



(2 - <$„,o) 



VJ4 



2i)e 



m 

/eB 



The case of freely streaming, untrapped neutrinos is obtained from the above 
equation by putting /i„ =0 and replacing \i e by {p, e + E w ). When the magnetic 
field is much weaker than the critical field for protons, only electrons are 
affected and the neutrino-nucleon scattering cross-section expression remain 
unchanged by the magnetic field. The numerical values however, are modified 
due to changed chemical composition. The cross-sections are given by 



Ovn{E v ) G|-coSg c 



V 



16tt 3 



{Cl + 3C 2 A )m* N 2 T 2 p e 



T 2 



exp- — T ' 



(57) 
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If neutrinos are not trapped, we get in the elastic limit 

^n(E„) G 2 f cos 2 , 2 

V ~ 16tt 3 ( v 3CA > mN 1 K 

The neutrino-electron scattering cross section is 



2G F C0S l ln 2 



which in the untrapped regime goes over to 



a ue 2G F cosg e frt2 

~ - 3^ ^ 



C\)&TEl 



(58) 



(59) 



(60) 



Non-Degenerate Matter : We now treat the nucleons to be non-relativistic 
non-degenerate such that [iijT <C — 1 and thus the Pauli-blocking factor 1 — 
fN{Ei) can be replaced by 1, the electrons are still considered degenerete and 
relativistic. The various cross-sections are given by 



v A {E v ,B) 
V 



G% cos 2 
2^~ 



-(C v + 3C 2 A )n n (E„ + Q) 



1 



1 _|_ e -/3(B„+Q-Me) 



1 



yJ[E v + Q) 2 -m 2 e - 2veB 



where tin is the nuclear density and Q — m n — m p . 

<J v n{E v ,B) G 2 F COS 2 g or i2\„ P 2 

y |Z \G V + 6C A )n N E v 

If the electrons too are considered non-degenerate we get 



(61) 



(62) 



1A 
V 

OvN 

V 



Gp cosg c 
G\ cos?, 



(giv + 3g 2 A) n NEl 



4k 



^(C v +3C 2 A )n N El 



(63) 



~N~ 



4Gp cos?, 



-T A E v {Cv +C A f 



(64) 



Quantising Magnetic Field : For quantizing magnetic field the square of the 
matrix elements can be evaluated in a straight forward way and we get 

\M\\ = %G 2 F cos 2 e c m* 2 ml{E i + Piz ) [{g v + g A ) 2 (E 1 + p u ) + 4g 2 A {E x - p u )] 
1 



exp 



2eB 



((Plx+P2x) + {P3x + Piy) ) 



(65) 
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\M\1 P 



16G 2 F C0Sg c C 2 A m 2 2 (pi ■ p 4 + 2p lz p 4z ) 
1 



exp 



2eB 



{{PAx +Plx) + {Piy ^Ply) ) 



(66) 



\M\L = 



16G 2 F cos 2 . 



(C'l + C 2 A )((E 1 E i + Pu p iz )(E 2 E 3 +p 2zP3z ) 
-{E lP4z + E 4 p lz )[E 2 p 3z + E 3 p 2z )) + 2C V C A ({E 1 E 4 + p lzP4z ) 
(E 2 p 3z + E 3 p 2z ) - Eip 4z + E 4 pi z )(E 2 E 3 + p 2z p 3z )) 
1 



exp 



2eB 



{(P4x -Plx) - (Piy-Ply) ) 



(67) 



The absorption cross-section is now given by 



V 



1 



d 3 p z 



eBL^/2 poo 



2E l L x J [2i:f2E 2 J_ eBL , 2 J_ QO (2i:f2E: 



dp 3y dp 3z 



eBL^/2 r oo 



dp 4y dp 4z 



eBLx/2J-oo (2TT) 2 2E 4 

f 2 (E 2 )[l-f 3 (E 3 )][l-U(E i )] 



(2ir) 3 6(P y )S(P z )S(E) \M\ 



(68) 



The corresponding scattering cross-sections are obtained by interchanging the 
particle 2 with 4. Performing the integrals with usual techniques we have 



<j a {E u B) eB 1 1 



V 



2Ei (2ir) s 



dE- 



dp 3z dp. 



•4z 



E 3 E 4 



6(E) \M\ 
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(69) 



To make further progress, we consider the cases of extreme-degeneracy 
and non-degeneracy separately. 

Degenerate Matter : For strongly degenerate matter, particles at the top of 
their respective fermi-seas alone contribute and in this approximation we get 
□ 

a A (E v ,B) G 2 F cos 2 6 c m* n m* p T 2 (n 2 + (^-) 2 ) 
-< /)• 
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exp 
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8\PF{n) +Pf(p) 



+exp[ — 



[ PF {n) + p F (v)] 2 - 4p|(e) 
2eB 
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Figure 1. Variation of Y p , Y n and m* with ns for B = 0, 5 X 10 5 and B K = 5 X 10 5 where 
in B K the effect of anomalous magnetic moment has been included. 



e[ \p F (n)+p F (u)} 2 -4p 2 F (e) 



(70) 



The case of freely streaming neutrinos is obtained by putting pf{v) = 
everywhere. 
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Non-Degenerate Matter: Similarly performing the integrals in the non degen- 
erate limit, the various cross sections are 



a A {E v ,B) G 2 f cos 2 9 c 
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eBcosg njv- 
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— ~ n e l hj v 

V TT 



{C v + C\){1 + cos 2 9) - ACvCacosO 



(72) 
(73) 

(74) 



We see that in the neutrino free case (Y Ve = 0), the neutrino absorption cross 
section i.e. the direct URCA process which is highly supressed for degener- 
ate matter in the absence of magnetic field, proceeds at all densities for the 
quantising magnetic field. 




Figure 2. Neutrino emmisivity as as function of Baryon density for B = and B = 10 5 , 
when the electrons and protons are completely polarized 



6 Results and Discussion 

We first calculate the composition of matter for arbitrary magnetic fields for 
both the neutrino free and neutrino trapped case over a wide range of density 
and temperature. We find that the effect of magnetic field is to raise the proton 
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Figure 3. Bulk viscosity as a function of Baryon density at T = 1 MeV for different values 
of magnetic field 



fraction and is pronounced at low densities. The effect of including anomalous 
magnetic moment also becomes significant at field strengths ~ 10 5 B^ which 
can be seen from Fig.l where we have plotted the proton fraction and effective 
nucleon mass as a function of density. We then calculate the effect of magnetic 
field on neutrino emissivity for the direct URCA process and find that the 
effect is not important at weak magnetic fields. However for the interesting 
case of quantising magnetic field capable of totally polarising the electrons 
and protons, we see from Fig. 2 that in this case the threshold for direct 
URCA process is evaded and the emissivity is enhanced by upto two orders 
of magnitude and develops anisotropy. 
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Figure 4. Variation of absorption mean free path with Baryon density. The first figure is 
for untrapped degenerate matter and the second figure for trapped degenerate matter. 
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7-ray burst events if they arise due to pulsations and the maximum ro- 
tation frequency are influnced by the bulk viscosity of neutron star matter 
through 1) damping of radial oscillations and 2) by influencing gravitational 
radiation reaction instability limiting maximum rotation rate. In Fig. 3 we 
plot bulk viscosity as a function of baryon density for different values of mag- 
netic field and find that the effect of magnetic field is very pronounced at low 
densities where only the lowest Landau level contributes and the direct URCA 
process is no longer inhibited. We also find that the viscosity decreases very 
rapidly with density. 

Neutrinos emitted during the cooling of protoneutron star have a mo- 
mentum ~ 100 times the momentum of pulsars and therefore, one percent 
anisotropy in neutrino emission would give a kick velocity consistent with ob- 
servation. This can occur in two different ways since the neutrino dispersion 
relations in magnetised medium are modified. Firstly due to neutrino scatter- 
ing reactions on polarised electrons and protons and secondly^ due to matter 
induced MSW v e <-> v T neutrino oscillations in the presence of magnetic field 
when v r sphere develops asymmetry and resonant surface gets distorted and 
the neutrinos escaping from different depths emerge with different energies. In 
Fig. 4 and 5 we have plotted the neutrino absorption and scattering mean free 
paths with density for different values of temperature and magnetic field for 
both the free streaming and trapped regimes. In Fig. 6 we show the variation 
of absorption mean free path with density and observe substantial decrease 
in addition to developing anisotropy with magnetic field which is particularly 
pronounced at low densities. 

To conclude, we find that the effect of quantising magnetic field and of the 
inclusion of anomalous magnetic moment of nucleons is to increase the proton 
fraction and lower the effective nucleon mass so as to evade the threshold for 
the direct URCA processes to proceed. This results in the enhancement of 
neutrino emission and bulk viscosity. This also has the effect of substentially 
decreasing the neutrino absorption mean free path in addition to developing 
anisotropy. 
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Figure 5. Variation of neutrino scattering mean free path with baryon density the left panel 
is for untrapped matter and the right panel is for trapped matter 
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